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Abstract 

The present status of theoretical and experimental investigations of the decay 
rate of a positronium is considered. The increasing interest to this problem has 
been caused by the disagreement of the calculated value of r3(o — Ps) and the 
recent series of precise experiments. The necessity of new calculations on the basis 
of the quantum field methods in bound state theory is pointed out with taking into 
account the dependence of the interaction kernel on relative energies. 
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The quantum-electrodynamic systems, consisting of particle and anti-particle, have 
specific features. Apart from a scattering channel, the annihilation channel appears in 
this case. A positronium atom, which is a specimen of these systems, has no stability. 
The life time of a positronium (or the decay rate) is the subject of precise experimental 
and theoretical investigations. The charge parity of a positronium, C = (—1)^"'"'^ (L is the 
eigenvalue of an angular momentum operator, 5* is the eigenvalue of a total spin operator 
for a system under consideration), is a motion constant. Consequently, all its states are 
separated into the charge - even states {S = 1) and the charge - odd ones {S = —1). The 
positronium total spin is also conserved and the energy levels are classified as the singlet 
levels {S = 0, parapositronium) and as the triplet ones {S = 1, orthopositronium) . The S - 
state {L = 0) parapositronium has a positive parity and the S - state orthopositronium has 
a negative parity. As a consequence of conservation of a charge parity in electromagnetic 
interaction a parapositronium is disintegrated into the even number of photons and an 
orthopositronium is decayed into the odd one. 

At the present time, the essential disagreement is turned out between the theoretical 
and experimental values for the decay rate of an orthopositronium. The theoretical pre- 
dictions are 
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-10.266 ±0.011, 
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Taking into account the modern value of a, the fine structure constant, the result can be 
recalculated 0: 

7.038 31(5) (4) 
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The last experimental measurings are 
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The result of Ref. has 9.4 standard deviation from the predicted theoretical decay 
rate and the result of has 6.2 standard deviation. The coefficient ^3 = 1 in the 
0{a^) term can contribute 3.5 ■ 10~^/is~^ (or 5 ppm of T^) only. To take off the above 
disagreement the coefficient has to be equal to about ~ 250 ±40, what is very unlikely. 



"See Table I for the preceding experimental results. 
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indeed. However, this opportunity is pointed out in |^ just not to be rejected a priori 
Therefore, the calculation of coefficient is desirable enough now. 

For the first time the main contribution in the orthopositronium decay rate has been 
calculated in ilH : 



To{o - Ps) = -2Im{AE: 
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The corrections of the 0{a) order of the magnitude to this quantity had been calculating 
by numerical method |0, |, |], || at first, but later some of them has been figure out 
analytically in [1^ and |TO[-|T^ in the Feynman gauge. The corrections, which come 
from the diagrams with the self-energy and vertex insertions, have been calculated by 
Adkins 0, IH: 
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is standard expression of a dimensional regularization {2uj is a space dimension). The 
indices "IV" and "OV" designate the insertions in the internal photon-electron vertex 
and in the outer ones, correspondingly. The above results are co-ordinated with the 
Stroscio's result iTOl when 
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being added to the last one. This is necessary because of the different regularization 
procedures which have been used in [|lOl and in |Tl|, |12|, correspondingly. 

Recently, the calculations of these corrections have been finished in the Fried - 
Yennie gauge: 
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The contributions from the remained diagrams ( with a radiative insertion in a vertex of 
an internal photon; with two radiative photons spanned; the diagram taking into account 
boundary effects and the annihilation diagram, see Fig. I in [5b] ), have been calculated 
numerically. Totally, the 0{a) corrections are joined to give 
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Then we have^ : 

F*^^g(o - Ps) = 7.038 236(10) /is-\ (19) 

The above result is the most precise theoretical result available at the present moment. 

To solve the existing disagreement between theory and experiment, the 5 - photon 
mode of o — Ps decay and the 4 - photon mode of p — Ps decay have been under consid- 
], p!5[]p[ The following theoretical evaluations have been obtained: 
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0.018 9(ll)a2r^, 
0.013 89(6)ma^. 

They are in agreement with each other and with the results of the previous papers [|IB[0 : 

[16]. o 



T^^ \p-Ps) = 0.013 52 ma^ = 11.57-10 
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(24) 

^The uncalculated yet 0{a^) corrections are not accounted here. 
As a consequence of a conservation of an angular momentum and an isotropic properties of coordinate 
space an orthopositronium has to decay into the odd number of photons and a parapositronium has to 
decay into the even one, as outUned before. 

'^The result is not correct, four times less than the above cited results. The explanation of this 
was given in m5 \ . 
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In the connection with the present situation with respect to the decay rate of o — Ps 
investigations of alternative decay modes for this system (e.g. o — Ps — » 7 + a, a is an 



axion, a pseudo-scalar particle with mass ma < 2me) are of present interest [T3]-[EH]. In 



the article |2T| the following experimental limits of the branching of decay width have 



been obtained: 

Er = ^- — < 5 • 10"^ - 1 ■ lO"*' 30ppm , 25 

r(o — Fs ^ 37) 

provided that rria is in the range 100 - 900 keV. In the case of the axion mass less than 
100 keV (which is implied by the Samuel's hypothesis [2^; according to the cited paper 



TJia < 5.7 keV, Qae+e- ~ 2 • 10 ^) the limits of Br are the following ones |Q : 



Br = 7.6 -10-^ if ~ 100 fceV, (26) 
Br = 6.4 ■ 10-^ if m„ < 30 keV. (27) 

These limits are about 2 orders less than the value which is necessary to remove the 
disagreement. 

Finally, a decay o — Ps — * nothing (that is into weak-interacting non-detected parti- 
cles) Q has been investigated in The obtained result 

r(o — Ps —>■ nothing) 4 



r(o -Ps^ 37) 



< 5.8 ■ 10"^ (350ppm) (2J 



expels the opportunity that this decay mode is an origin of disajustment between theory 
and experiment. 

The decay of o — Ps into two photons, which breaks the CP - invariance, as else 
mentioned in [26, was experimentally rejected in |28] 0. 



Let us mention, the contribution of weak interaction has been studied in ||30|1. However 



because of the factor m^/M^ ~ Gp ■ frte — 3-10 it cannot influence final results. In 
the cited article the weak decay modes are estimated as 

r(p-Ps->37) r(o-Ps^47) 2 2 in-27 

r(p-p.^27) " r(o-p.^37) " " ^0 ' (30) 

where Gi? is the Fermi constant for weak interaction, 

gv = l- ^sin^Qw ^ 0.08, (31) 

*The proposed values do not influence the agreement of theoretical and experimental results of an 
anomalous magnetic moment (AMM) of an electron. 

•^Analogously to Glashow's hypothesis of the decay into invisible "mirror" particle p^ . 

f The physics ground of these speculations is a possible existence of an unisotropic vector field with 
non-zero vacuum expectation , with whom an electron and a positron could be interacting 

C = giPOaf3^A"n'^, (29) 

C is the interaction Lagrangian. 
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Qw is the Weinberg angle. The current experimental limits are 

T{p-Ps 37) 



T(o 



Ps 
Ps 



27) 
47) 



r(o-Ps ^37) 



< 2.8 -10-^ 
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In the Table I all experimental results for the o — Ps decay rate , known to us, are 
presented 0. 

Table I. The experimental results for the o — Ps decay rate. 



Year 



1968 
1973 
1973 
1976 
1976 
1978 
1978 
1978 
1978 
1982 
1987 
1987 
1989 
1990 



Reference 



ii 



m 

1 38 



39| 
40| 
111 



m 



m 



111 



3B| 



m 



13(0 - Ps), flS ^ 



7.262(15) 

7.262(15) 

7.275(15) 

7.104(6) 

7.09(2) 

7.056(7) 

7.045(6) 

7.050(13) 

7.122(12) 

7.051(5) 

7.031(7) 

7.0516(13) 

7.0514(14) 

7.0482(16) 



Error, ppm 



2070 
2070 
2060 
840 
2820 
990 
850 
1840 
1680 
710 
1000 
180 
200 
230 



Technique 



gas 
gas 
gas 

powder 5*^02 
vacuum 

gas 

gas 
vacuum 
vacuum 

gas 
vacuum 

gas 

gas 
vacuum 



Regarding the results for the decay rate of a parapositronium, the situation was highly 
favorable until the last years. The theoretical value, which was found out else in the fifties 



48, 49 , is equal to 



T^ip - Ps) = -2Im{AE. 



27 J 



1 a^mc^ 

2 h 



TT^ 4 



7.9852 ns-\ 



(34) 



The above value, confirmed in [5C, 51], coincides with the direct experimental result up 
to the good accuracy: 

rJlgj(p-Ps) = 7.994 ± 0.011 (35) 
The experimental values of the parapositronium decay rate are showed in the Table 11 []. 



''The results of the papers [g3[ Q and Q can be accounted as rough estimations only. 

*The branching of the decay rates of a para- and an orthopositronium p^|^~p^j had been measuring 
in the experiments of 1952 and 1954 only. The presented results are recalculated by means of the first 
direct experimental value T3{o — Ps) = 7.262(15) /is~^ Bq]. 
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Table II. The experimental results for the p — Ps decay rate. 



Year 


Reference 


T2{p — Ps), ns ^ 


Error, % 


Technique 


1952 




52 




7.63(1.02) 


13 


gas 


1954 




53 




9.45(1.41) 


15 


gas 


1970 




54 




7.99(11) 


1.38 


gas 


1982 




45 




7.994(11) 


0.14 


gas 



In the articles [y, |5T| it was pointed out that it is necessary to add the logarithmic 
corrections on a to the Harris and Brown's result. In the article [13a] these corrections to 
the T^i^o — Ps) and T2{p — Ps) have been calculated again, with the result of the decay 
rate of a parapositronium differing with the one found out before p]] : 

T^^\p- Ps,aHna) = ■ 2aHna-\ (36) 

V^^^{p-Ps,aHna) = .-aHna-\ (37) 

2 3 

Finally, the quite unexpected (and undesirable) result, presented in the Remiddi's 
(and collaborators) talk |^ should be mentioned. The calculations, carried out by the 
authors of cited paper, lead to the additional contribution: 



rPl(p - Ps, alna) = ^(-)2/n«, (38) 



which is explained by authors to appear as a result of taking into account the dependence 
of an interaction kernel on the relative momenta.^ 

The above-mentioned shows us at the necessity of a continuation of calculations of 
the decay rates of an orthopositronium as well as a parapositronium employing more 



■'Let us mark, the additional contributions, which are similar to the ones of Ref. ||5q| , appeared 
in the calculations of the hyperfine splitting (HFS) of the ground state of two-fermion system by the 
quasipotential metho d |56| [5^ when taking into account the dependence of the interaction kernel on the 
relative energies ^9|. For example, the 0{a^lnQ) correction to the HFS, obtained from one-photon 
diagram, is equal to 

AE^/^iaHna) = Ep-^^C^ + ^ + 2)lna-\ (39) 
mim2 m2 mi 

when using the version of the quasipotential approach based on the amplitude t = 

(G|^)^^GoTGo (Gq)^^ ■ Here G+ is the two-time Green's function projected onto the positive-energy 
states, mi and m2 are the masses of the constituent particles, Ep is the Fermi energy. The index " fr^' 
designates that the diagram of one-transversal-photon exchange is under consideration. 

Using the kernel constructed from the on-shell physical amplitude we can find out that the last term 
in the brackets disappears. However, the total result to the HFS accounting all diagrams is the same in 
both of versions of quasipotential approach. 

We can also come across with the similar situation in calculations of the anomalous corrections of 
0{alna) order, which do not appear in the method on the mass shell. However, the version, based on the 
two-time Green's function, shows the contraction of these anomalously large terms, taking into account 
two-photon exchange diagrams [^9| . 
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accurate relativistic methods, one of which is the quasipotential approach in quantum 
field theory |5^ giving the opportunity taking into account of binding effects. 
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